Resonant trapping in the transport of a matter-wave soliton through a quantum well 
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We theoretically investigate the scattering of bright solitons in a Bose-Einstein condensate on 
narrow attractive potential wells. Reflection, transmission and trapping of an incident soliton are 
predicted to occur with remarkably abrupt transitions upon varying the potential depth. Numerical 
simulations of the nonlinear Schrodinger equation are complemented by a variational collective 
coordinate approach. The mechanism for nonlinear trapping is found to rely both on resonant 
interaction between the soliton and bound states in the potential well as well as radiation of small 
amplitude waves. These results suggest that solitons can be used to probe bound states that are 
not accessible through scattering with single atoms. 
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I. INTRODUCTION 

A classical particle incident on a potential barrier can 
either reflect or pass over it, depending on its kinetic en- 
ergy in relation to the height of the barrier. Therefore 
it always passes a negative barrier, e.g. a hole. This 
is not the case any more for quantum mechanical parti- 
cles or matter waves [l|. Solving the linear Schrodinger 
equation shows that partial reflection can be expected 
from attractive potentials such as quantum wells. In re- 
cent experiments by Pasquini et al. 0, Q Bose-Einstein 
condensates (BECs) have been found to reflect from a 
surface in spite of mostly attractive atom-surface interac- 
tions. In addition to the wave nature of ultra-cold atoms, 
BECs also experience nonlinear mean-field interactions, 
which present a complication in the experiments 0, [1[ . 
On the other hand, the nonlinear interactions are po- 
tentially useful when they result in effects that are not 
obtainable with linear matter waves. 

Here we focus on BECs with attractive interactions 
confined to a quasi-one-dimensional waveguide that can 
form self- localised wave packets known as bright solitons. 
In the limit of tight transverse confinement, the mean- 
field theory of BECs reduces to the nonlinear Schrodinger 
equation, an integrable soliton equation Q. A well- 
known property of solitons is that they behave in many 
respects as classical particles do, e.g. with respect to their 
collisional properties or their motion in a slowly-varying 
external potential [?J • As solitons have both particle and 
wave properties, they may experience "quantum" reflec- 
tion from an attractive potential well and yet maintain 
their particle- like integrity to a large degree @. In ad- 
dition to such nonlinear wave effects, macroscopic quan- 
tum tunneling and fragmentation have recently been dis- 
cussed UGH. 

On the other hand there is a possibility that a scat- 
tering soliton, or part of it, may become trapped in a 
potential well. Such an effect would not be possible for 
either a classical or a quantum mechanical particle in the 
absence of dissipation. In this paper we study the effect 
of trapping of solitons in attractive potential wells. We 



argue that a resonant population transfer between the 
soliton and non-linear bound states of the potential well, 
first suggested by Goodman, Holmes, and Weinstein [ll|, 
are vital for this effect. Here, we explore the details and 
the consequences of this mechanism. We find that it is 
possible to deduce the energy of the most weakly bound 
state of an unknown defect by scattering solitons with 
known parameters and by recording the trapped particle 
number. This could potentially lead to real-world appli- 
cations of nonlinear-wave scattering. Although the cur- 
rent work addresses matter- wave solitons in particular, 
our findings are equally applicable to nonlinear optics or 
other nonlinear wave problems governed by the nonlinear 
Schrodinger equation. 

Recent progress in experiments has made it possible 
to create matter- wave solitons and to explore their prop- 
erties [III, [TH, [HI- This, besides the importance of soli- 
ton physics in other areas of physics [la, [H, H3; Gil has 
motivated a variety of authors to investigate the scat- 
tering of solitons on different kinds of potentials, like 
barriers or impurities [IMlJSJ ES EH El El, wells 
IJHLElESlli], steps TUlill and a potential ramp 
301 ] - Here, we focus on the quantum well, as which we 
understand an attractive potential well with well distin- 
guished single-particle energy levels representing linear 
bound states. One of us has previously investigated the 
enhancement of quantum reflection by nonlinear interac- 
tions in solitons and the abrupt transition to transmis- 
sion in the scattering on an attractive defect potential 
0. The aim of the current work is to extend the pre- 
vious work to include trapping phenomena. Varying the 
strength of the attractive defect and scattering a slow 
soliton, we identify regimes dominated by transmission, 
reflection, trapping, and a combination of trapping and 
reflection. We discuss a trapping mechanism by a reso- 
nant transfer to (quasi-) bound states within a quantum 
well. Similar resonant effects have been investigated in 
the transport of repulsive condensates through a double 
barrier potential [31] . 

In Sec. Ullof this paper we introduce the theoretical ap- 
proach employed for scattering matter-wave bright soli- 
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tons in a tight waveguide trap. Numerical simulations 
of trapping phenomena are presented in Sec. IIIII before 
discussing resonant mechanisms for trapping and trans- 
mission. A collective coordinate approach based on a 
variational model is discussed in Sec. IIVI and compared 
with the simulation results. Sec. [V] then discusses the 
trapping mechanism in more detail. The final section Sec. 
IVTl discusses how energy levels of defects can be probed 
via scattering of solitons. 



II. THEORETICAL MODEL 
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We consider an attractively interacting BEC in a 
waveguide-like trap with tight harmonic confinement in 
two dimensions but weak or no confinement in the re- 
maining (z) dimension. In such a situation, bright soli- 
tons constituting localized (bound) BEC wave packets 
are metastable and collapse ensues beyond a critical par- 
ticle number [32|. We assume that the soliton size re- 
mains sub-critical and that the linear density n(z, r) at 
any time remains well below the threshold for transverse 
collapse — na s <C 11.7/(87r) [H, [34[, where a s is the 
s-wave scattering length. In this case, the soliton dy- 
namics may be modeled with the one-dimensional Gross- 
Pitaevskii (GP) equation 



n 9 a> 
iri—cp = 

or 



4>, (i) 



where the GP wave function \4>{z, r)| = \J n(z, r) is nor- 
malized to the number of atoms J |0| 2 <iz = N and the 
ID interaction constant is g\u « 2hajj_a s [35L [36| under 
the influence of a transverse harmonic trapping potential 
with frequency uj±. The quasi-one-dimensional approxi- 
mation {1} is expected to break down for solitons of close 
to critical size, where the dynamics may become inher- 
ently three-dimensional [371 ] ■ Small quantitative correc- 
tions to Eq. |T]) due to the finite transverse extent of the 
solitons (33, [39| are not expected to significantly alter the 
results reported below and are thus neglected. 

Choosing an energy scale E > and a density scale h 
we can rewrite Eq. (fl]) in dimensionless form by introduc- 
ing t = r/t, x = z/x, i){x) = (f>(z)/\/fi, where t = h/E 
and E = ti 2 / (mx 2 ): 



9 << \ 



l d 2 

2dx^ 



+ g\^(x,t)\ 2 + V(x) 



ip(x,t). 



(2) 

where we have introduced V{x) = V\d{x)/E and the 
dimensionless coupling constant g = giDn/E. In the 
following we assume g < since we consider attractive 
BECs that support bright solitons. At this point the 
energy E and density scale n remain arbitrary and can be 
chosen to suit experimental parameters. We will discuss 
specific choices below. 

For a vanishing potential V(x) Eq. ^ has the soliton 



FIG. 1: A soliton is being scattered on a rectangular well with 
width 2a and depth Vb. 



solution 



ip(x) = Asech (A A 



x Q — vt)) exp(i6>(a;, £)), (3) 



where v is the velocity of the soliton measured in units 
of x/t and xo is the dimensionless position at t = 0. The 
solution is normalized according to 



Af s = / dx\i;(x,t)\ 2 = 



2A 



(4) 



where A is a dimensionless amplitude and J\f s is related 
to the particle number by N — N s nx. The phase is given 
by 9{x, i) — vx — ujt and uj — t> 2 /2 + /i is the dimensionless 
frequency. Here, /i = gA 2 /2 is the (negative) chemical 
potential of a stationary soliton measured in units of E . 
The soliton width in units of x is given by l s = l/(AJ—g) 
0. 

For the numerical simulations below we will commonly 
choose g — —1, which relates the energy and density 
scales by E = gmn. The further choice of A = 1 
fully determines the energy scale E = mN 2 g 2 D /(Ah 2 ) = 
N 2 <J 2 _a 2 m and the density scale n = N/(2x), where the 
unit length becomes x = —h/(Nujj_a s m) and the dimen- 
sionless soliton length becomes l s — 1. 

Typical experimental values for 7 Li BEC [13, EH are 
luj_ ps 2tt x 710Hz, a s -0.2nm and N w 6 x 10 3 . This 
yields a length scale of x = 1.7 um, which is consistent 
with experimental observations [12]. For these param- 
eters the time unit is t w 0.3ms. The velocity scale is 
consequently x/t « 5.7/xm/ms. 

For the form of the external potential we consider a 
rectangular well defined as 



V(x) = 







for \x\ > a 
for |a:| < a. 



(5) 



as shown in Fig. [T] We are specifically interested in the 
case where the width of the well 2a is comparable to the 
soliton width / s and therefore we choose 2a = 1 for the 
numerical studies in this work unless noted otherwise. 
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Lee and Brand [8| have already investigated in detail 
the enhanced reflection of solitons for the special case of 
a Rosen- Morse potential — Vbsech 2 (ax) at low velocities. 
There they found a step-like behavior of the reflection 
and transmission probabilities, which is due to the non- 
linearity. Here, we aim to extend this work to include 
resonant trapping effects. 

We solved Eq. @ numerically via the Crank-Nicholson 
method using a standard finite difference discretization 
of the spatial derivatives • The algorithms were im- 
plemented in standard C/C++ and Octave [4l|. We per- 
formed the simulations in a box with hard wall bound- 
aries. The box length was set to Ibox = 80/ s unless stated 
otherwise. Furthermore we used N g = 2001 grid points 
and a fixed time step of At = O.Of. The convergence 
of our calculations with respect to these quantities was 
monitored carefully. Reflection from the boundaries was 
avoided by appropriate timing the the simulation. We 
also used complex absorbing potentials at the bound- 
aries for verifying that reflection effects remained below 
a quantifiable threshold. 




Potential Depth V 



FIG. 2: Reflection R, Transmission T and Trapping L (L for 
'L'ocalized) as a function of Vo and a fixed velocity v — 0.3 
by solving Eq. [2] The lower picture zooms into one of the 
structures shown in the upper picture. The well width 2a = 1 
is kept constant and we use A — 1 and g = — 1. The same 
parameters are used throughout the paper unless explicitly 
mentioned otherwise. 



III. PHENOMENOLOGY OF SOLITON 
SCATTERING ON A QUANTUM WELL 

In this section we present results from numerical solu- 
tions of Eq. ([2]) corresponding to a soliton approaching 
the well of Eq. ([5]). In the initial setup the soliton © is 
being placed at position xq = —12 left of the quantum 
well moving with the velocity v > towards it (see Fig. 

As physical observables we introduce the reflected (R) , 
trapped/localized (L), and transmitted (T) fraction of 
the soliton, which are calculated at a time well after the 
initial impact of the soliton on the well (e.g. t — 166): 



R 
L 
T 



1 

KJ- 



dx\ip(x, i)| 5 



oo 
8o 



1 f 

1 f°° 

w / dx\iP(x,t)\ 2 



(6) 



with R + L + T = 1. Figure [2] shows these quantities 
as a function of the depth of the well for a fixed initial 
velocity ^initial = 0.3. We consciously study the case 
of small velocity where v? nitial /2 <§; |/x| = 0.5. For the 
parameters of Ref. [l3[ (see also Sec. HI]) , this velocity 
amounts to « 1.7mm s . 

The upper panel of Fig. [2] shows several structures of 
similar form on a background of almost complete reflec- 
tion. We thus call these structures reflection-trapping 
(RL) windows. In the following, we focus our discussion 
mostly on the second one as shown in the lower panel. 
For a certain range of Vo the soliton reflects completely on 
the well. But by increasing the depth of the quantum well 
R suddenly drops to zero while the transmitted fraction 



jumps to an absolute maximum. Further increase gives a 
sudden drop of T to almost zero and most of the soliton 
is being trapped inside the quantum well. Then the trap- 
ping component L starts to decrease while the reflected 
part increases. At least some of the reflected and trans- 
mitted amplitude in this part of Fig. [5] can be attributed 
to radiation, i.e. small amplitude waves. This becomes 
apparent in Fig. [3J where snapshots of the density of the 
time-dependent wave function are shown. We will dis- 
cuss the role that radiation plays in enabling trapping 
by carrying away kinetic energy in Sec. [VJ At slightly 
larger Vq we observe the co-existence of a reflected soli- 
ton with a trapped component together with radiation in 
the transmission channel becoming very small again. 

Figure [3] reveals another remarkable feature: The con- 
densate density has a single node localized close to the 
center of the well. Our simulations show that the num- 
ber of nodes located in the well is a characteristic of each 
RL window. Indeed we find that RL windows appear 
around a critical well depth, where a linear bound state 
with the appropriate number of nodes is formed. In the 
first RL window, the density reveals no node, the second 
one shows one, the third one shows two nodes and so on. 
The density of the soliton while located over the well is 
similar to the density functions for bound states of the 
Schrodinger equation in a quantum well. A more detailed 
analysis of the relation of the RL window to linear res- 
onances and nonlinear bound states of the well is given 
further on in this section. 

The time-dynamics of the soliton are summarized in 
the density plot in Fig. |4j The pictures show the four 
different scenarios of full reflection, full transmission, full 
trapping and partial trapping. On the lower left picture 
the density sloshes around the center but a closer look 
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FIG. 3: The upper picture shows a snapshot of the condensate 
density at t = 77 for Vo = 5.2 where trapping is maximized 
in the second RL window. One can see the trapped mode 
in the first excited bound state (see text) and the radiation 
which stabilizes the trapped soliton. The lower picture shows 
partial trapping with a reflected soliton at t=65. 



brings up that the radiation reduces the amplitude of this 
oscillation and therefore stabilizes the trapped soliton. 
Furthermore, the position of the dip in density remains 
almost stationary, varying by not more than 5% of the 
potential width. The reason for this is that for our choice 
of parameters the energy differences between the bound 
states in the well are large compared to any energy scale 
of the incoming soliton. Hence only one of these states 
can be populated, in case of Fig. [4] it is the first excited 
state. 

We now discuss the relation of the trapping phe- 
nomenon to (stationary) nonlinear bound states of the 
well. Figure [5] compares two different observables. The 
first one is the trapped component L from the time- 
dependent simulations. The other one gives the relative 
number of particles Nl^^Vq) in an eigenstate of the 
time-independent Gross-Pitaevskii equation for a fixed 
chemical potential /z, that is set to the same value as the 
chemical potential of the free soliton (initial = —0.5) in 
the time-dependent simulations. Specifically, -/Vt jre ;(Vo) 
is given by 



N, 



L.rel 



(Vo) = 



(7) 



where A/e(Vo, fx) is the normalization constant (HJ) of the 
single-node stationary solution of Eq. ((2]) with the chem- 
ical potential /x while Afs(fJ-) is the normalization of a 
free soliton with the same chemical potential. In the nu- 
merical procedure Vo is changed iteratively to keep the 
chemical potential at the desired value. The results for 
NL,rei(Vo) can then be compared with the relative num- 
ber of trapped atoms L we got from the time-dependent 
simulations. Even for different parameters the agreement 
between both graphs is very good. These findings indi- 
cate that trapping is a resonant phenomenon with the 
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FIG. 4: Time and spatial dependence of the condensate den- 
sity \ip(x,t)\ 2 in gray scale (normalized to a maximum ampli- 
tude of 1) for four different Vo as in Fig. [2]but with lt ox ~ 40 
The case of full reflection is shown in panel a) while the one 
for full transmission is given in b). Furthermore c) presents a 
fully trapped soliton while in d) the case of partial trapping 
and reflection is shown (the additional reflection towards the 
end comes from the hard wall boundary conditions). 



chemical potential being the parameter of primary rele- 
vance. 

Another feature in Fig. [2] are the resonant transmis- 
sion bands. They are closely related to the above-barrier 
transmission resonances in the linear Schrodinger equa- 
tion, which is found from Eq. [T]for g = 0. There one can 
find the analytical solution for the transmission [42j 



Tun(Vu) = 



1 



V 2 
1 o 



■ sin 



; (2a v / f 2 + 2V ) 



v 2 {v 2 + 2Vo) 

In Fig. [5] we compare the transmission for g — with the 
case of solitons at g = — 1 at different velocities v (see 
Fig. 0. For very high velocities both curves approach 
each other. This is easily explained by the fact that the 
kinetic part in Eq. @ becomes much larger than the 
nonlinear term and therefore dominates the transmission 
spectrum. Thus decreasing v increases smoothly the ef- 
fects of the nonlinearity, in particular the formation of 
resonant transmission windows instead of transmission 
resonance lines. But their positions remain the same, 
which means that the nonlinearity just affects the shape 
of the transmission lines. We conclude that the basic 
mechanism of above-well shape resonances known from 
the linear Schrodinger equation remains valid for solitons. 



IV. VARIATIONAL ANSATZ 

Goodman et al. [ll[ studied soliton-defect interactions 
by simple two-mode models featuring a mobile soliton 
and a localized (trapped) mode. Here we extend this 
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approach by including breathing of the trapped mode. 

We approximate the well by an attractive delta poten- 
tial, defined as 



V{x) = -S(x)Vo 



(9) 



V.. 



V.. 



with Vq > 0. With this simplification there is exactly one 
linear bound state for all potential depths. Therefore we 
use an ansatz that splits the total wave function 



FIG. 5: L(Vq) (solid line) from the time-independent solu- 
tions and NL,rel(Vo) (dashed line) from the time-independent 
calculations for a fixed chemical potential y, = — |. Both 
quantities show similar behavior, even for different potential 
widths and different bound states. 



into a free soliton 
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FIG. 6: Comparison of T for solitons (g = —1, solid line) 
with the analytical solution Tu n (Vo) [Eq. jHJ] for linear waves 
(g = 0, dashed line). From top to bottom we increased Uinitiai 
of the incoming soliton. For increasing velocities the kinetic 
term in the Gross-Pitaevskii equation becomes dominant and 
therefore both curves approach each other. 



ip s = A s sech(A s x — Q s )e l "e 



i$ s „iV B x 



and a trapped part 



ipt = ^tsech(a;/a t )e l * t e i<Ttlogcosh(2:/at) 



(11) 



(12) 



that models a nonlinear mode that is localized at the well. 
Here we introduced a particular form of chirping term 
logcosh(x/a t ), which is capable of describing breathing 
modes. This can be used as a substitute for radiation 
effects which should allow the soliton to be trapped as it 
can transfer kinetic energy into another form of excita- 
tion. The choice for this particular form of the chirping 
term is consistent with Ref. (24|. This leads to our sys- 
tem's Lagrangian given by 



L = 



dx 



+ -\^-V(x)\ 
1 A? 



-2Aia t $t - 2A\a t G t {2 - log(4)) + 2A 1 a t a t - --*■(! + of) + -A\a t 



3 a t 



-2A S $ S - 2V S Q S 
+V {A 2 t + A 2 s sech 2 {Q s ) + 2A t ,4 s sech(Q s ) cos($ s - * t )} 



(13) 
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To obtain the equations of motion one has to solve the for = A s , $ s , Q s , V a , A t , 3>t, a*, ot. 
Euler-Lagrange equations 

dt{wJ = W i (14) Thlsleadsto 
i 



— A s = V A s A t sech(Q s )sm{§ s -<S> t ) 

= \{A 2 a -Vf) + V [A s sech 2 (g s )+^sech(Q s )cos($ s -$ t )] 

~t:Qs — A S V S 
at 

^-V s = -V [^sech 2 (g s )tanh(g s ) + A s A 4 sech(Q s )tanh(Q s ) cos($ s - 

^o* = ^ + (2-log(4))^A s sech(g s )sin($ s -$ t ) 

at 6at A t 

^crt = i(l + o 2 ) - ^ - ^ + ^-^8ech(Q.) [2<7 t sin($ s - <f> f ) - cos($ s - $ t )] 

j t $ t = ~& t {2 - log(4)) + ^ - g^j (1 + + §4? + 2^ t^t + v4 s sech(g s ) sin(<I> s - $ t )] . (15) 
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time t 

FIG. 7: The phase difference A"I> for different potential depths 
and for the time when the incoming soliton reaches the quan- 
tum well at t ~ 30 from solutions of Eqs. (115[) . The solid line 
and the dotted line show clearly a phase difference close to 
where transmission respectively trapping occurs. An example 
for reflection with a turning point at t ~ 29 is given by the 
dashed line where A$ ~ ±7r (for initial conditions see text). 



Two of these eight dynamical variables can be elimi- 
nated due to conservation laws. The amplitude A t can 
be found from the normalization of the wave function 



Wo = 2A S + 2A\a u (16) 



since N§ is a constant of the motion. We obtain 

a, . MZE. (17, 

V 2 a t 

Furthermore, it is not necessary to calculate <£>t and $ s 
separately as the only interesting and physically impor- 
tant property is the phase difference 

A$ = $ t -$ s . (18) 

We note that the total phase <& t + <E> S has a trivial time 
dependence as it is canonically conjugate to the total 
energy (Hamiltonian), which is a further constant of the 
motion. This way the problem left to solve consists of six 
coupled first-order ordinary differential equations. 

This ansatz allows us to calculate the time dynamics 
of a soliton without solving the Gross-Pitaevskii equation 
directly. But, of course, this is still a very simple approx- 
imation and thus the results are not expected to be as 
accurate as the GP results. However they can give further 
insight in the mechanism involved. Like in the previous 
section the initial velocity is set to initial = 0.3 and the 
initial position of the soliton is Q s = —10. We choose 
physically reasonable, small initial values for the param- 
eters of the trapped mode (A t = 10~ 4 and at = 10~ 2 for 
the simulations) in order to avoid numerical divergences. 
Furthermore we use A$ = at t = but we find that 
the results do not depend on this initial choice. 

Figure [7] shows the phase difference A<£> for different Vb 
during the scattering process. For Vq = 0.2 the soliton 
is being transmitted and we find A<!> <C 1. Increasing 
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FIG. 8: Results for soliton scattering on a quantum well from 
the collective coordinate Eqs. (|15p . From top to bottom: Re- 
flection R and transmission T vs Vq, the total trapping Ltotal, 
fraction Lt in the trapping mode, fraction L s in the soliton 
mode. 
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FIG. 9: (Color online) Condensate density as a function of 
time by solving Eqs. (|15[l . analogous to Fig. [4] For Vb = 
0.2 panel a) shows the typical situation a fully transmitted 
soliton. Partial trapping at Vb = 0.7 is shown in b). Panels 
c) and d) present the cases for Vb = 0.38 and Vb = 0.78. 
The soliton is being trapped for c) due to the population of 
the trapped mode and sloshes around the well for d) while 
populating the soliton mode only. The last plot for Vo = 1.75 
shows the whole soliton being completely reflected. 



the potential depth to Vb = 0.5 results in (partial) trap- 
ping and A$ -C 1 while for a even larger Vb = 2 there 
is full reflection (A$ sw ±7r). From Eq. ((15]) we can see 
that the time dependence of the velocity V s highly de- 
pends on cos(A$). There we find that for a small phase 
difference the velocity does not change sign for all times 
and therefore the soliton transmits through the well or 
gets trapped. In the other case of A<E> « ±7r the soliton 
can reflect from the well as the sign of the velocity can 
change. The difference between trapping and transmis- 
sion however lies in the potential depth that determines 
how fast the trapped mode can be populated, i.e. large 
values for Vb result in a faster population as can been 
seen in Eq. (JT5J for A s and dt- Therefore we can find a 
band between the reflection and the transmission regime 
where trapping can occur. 

It is a well-know feature in collisions between bright 
solitons that a 7r-phase difference induces repulsion [43 1. 
This mechanism for reflection was discussed in Ref. [8|. 
In particular it avoids trapping. Conversely, a resonant 
process with small phase difference is responsible for the 
population of the trapped mode. This is consistent with 
the findings of the previous section where trapping was 
described as a resonant process. 

To complete the comparison with the previous section, 
Fig. [5] shows the reflection i?, trapping L and transmis- 
sion T as a function of the potential depth Vo- We see 
similar features as in FigO For very small Vb the soliton 
is being transmitted almost completely (Fig. (9^) while 
for large Vb full reflection (Fig. (U) can be observed. Be- 
tween both of these regimes we find a more complicated 
and interesting behavior. There, (almost) all the time 
one observes partial trapping of the soliton at the end of 
the simulation. Furthermore, we find two forms of trap- 
ping. The first case is the normal one. There we see 
that the trapping mode is being populated by the incom- 
ing soliton. The other fraction that remains in the soliton 
mode is moving either to positive or negative infinity (see 
Fig. [5]d and (5fc). In addition, another kind of trapping 
can be observed. In this situation the soliton mode os- 
cillates around the delta potential (Fig. (5]i). According 
to the numerical simulations this is the only event when 
full trapping occurs. 

We conclude that the basic ideas from the previous sec- 
tion are still valid: For small Vb there is full transmission, 
then (partial) trapping and for very large Vb the soliton 
reflects completely. 



V. THE TRAPPING PROCESS 

In order to study the role of energy conservation and 
radiation in the trapping process, we consider the energy 
functional 



E[ip(x)} = / dx 



v>Mz)l 2 + fhMz)l 4 

(19)" 
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FIG. 10: Transmission T rc i(K>) = T(Vb)/T(4.85) and ve- 
locity contribution to the kinetic energy E^ in rol (Vb) = 
-Ekin(Vb)/£km(4.85) relative to the values at Vo = 4.85 where 
transmission is maximal. As the transmitted fraction of the 
incoming soliton decreases for deeper wells, so does its veloc- 
ity Vt- 



FIG. 11: Time dynamics of different kinetic and interaction 
energies for the transmitted part of the soliton for Vo = 4.9 
(solid line) and Vo = 5.1 (dashed line) well after the collision 
(t « 45)). The soliton undergoes breathing oscillations after 
the scattering process that shown up as oscillations in the en- 
ergy. For deeper wells E% in becomes smaller, i.e. the velocity 
of the transmitted soliton v t is smaller. 



We split this energy into different energy terms 
These are defined as 



E kin + E kin + E int ■ 



(20) 



E 



kin 



dx 



dx 



d\i>{x)\ 



dx 



\ip{x)\ — ex${i6{x,t)) 



— I dx 



-Mx)\ A 



(21) 



with i/j(x) = \ip(x)\ex-p(i8(x,t)). The first term gives 
the contribution to the kinetic energy from the density 
variations while the second term represents a contribu- 
tion from the phase gradient, which is connected to the 
superfiuid velocity [6|]. Ei n t is the interaction energy. 
Specifically for the soliton solution Eq. we find 



E. 



kin 



l 



A 



-V , Ei n t — — — 



and for the total energy 



E[i> s ] = -^J-A 3 



A 



3 V-9 
In particular, we find the universal ratio 

Eint 



kv. 



= -2 



-gA 3 

(22) 
(23) 

(24) 



for the soliton solution. 
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FIG. 12: Logarithmic density plot of the condensate for dif- 
ferent times at Vo = 4.9. A reflected component (radiation) 
is clearly visible. 



We now show that radiation loss during a scattering 
event leads to a decreased velocity due to energy conser- 
vation. We consider a soliton @ with initial velocity Vi 
that, during a collision event, suffers a small loss in ampli- 
tude due to radiation (small amplitude waves spreading 
away from the soliton) . The amplitude is reduced by the 
effect of radiation to A t = A — e with < e <C A The 
energy of the transmitted soliton traveling with velocity 
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FIG. 13: Logarithmic density plot ol the condensate for dif- 
ferent times at Vb = 5.1. Compared to Fig. [12] the reflected 
part clearly larger. 



FIG. 15: The plot shows the maximum amount of atoms L max 
that can be trapped around the second transmission reso- 
nance. L max drops significantly with increasing velocity as 
proposed by Eq. ([25)) . 




time t 

FIG. 14: Time dynamics of different kinetic and interaction 
energies for the reflected part of the soliton for Vb = 4.9 and 
Vo — 5.1. The ratio LEi a t/-Eki n <C 1, unlike what one expects 
from a soliton. As the interaction energy is small the reflected 
part is mainly radiation. This is only valid in the transmission 
window the reflected fraction is a soliton again for the case of 
trapping 



Vt is given by 



EtbP) = 



I / a \3 (A _ e ) 2 

-(A - e) d + y 




5v + 0{e z ) + 0(5v 2 ), 



(25) 



where the last line has been linearized in e and Sv = 
Vt — Vi. Identifying the term in square brackets as the 
energy of the initial soliton and assuming that radiation 
loss carries away a positive amount of energy (since the 
only negative contributions to energy could come from 
the nonlinear term, which is assumed to be small for 



radiation), we realize that the linear term in Eq. 
must be negative. This leads to 



vt < Vi 



2vi 



e < Vi 



(26) 



since e > 0. The slowing down of solitons after the colli- 
sion can be seen in Fig. 1101 which compares the velocity 
part of the kinetic energy Ef. in oc vf and the transmis- 
sion. 

The same parameters are used as in Fig. f2J where the 
transmission window was found between Vq — 4.85 and 
Vb = 5.2. Fig. ITD1 shows that a small change in the trans- 
mission (w 5%) results in a strong decrease of the trans- 
mitted soliton's velocity vt (~ 70%). Extrapolating Eq. 
([26]) beyond the regime of small e we find that for 



2± 

-gA 



(27) 



it predicts v t ~ 0, which allows the soliton be trapped in 
the well. 

For Vi > J -§-A we find that the right hand side of 

Eq. (|26| is always positive as e < A. Therefore we expect 
that trapping is reduced until it vanishes for very high 
velocities Vi when kinetic energy dominates over nonlin- 
ear energy contributions. Then the system becomes ap- 
proximately linear and can be approximated by a single- 
particle. 

In Figs [T0lfT4l we show results for the energy contribu- 
tions after the soliton-well collision. In particular, Fig. 
im gives energy contributions of the transmitted part of 
the soliton after the collision. There we can find that 
for both cases -^r 1 ~ — 2. The curves show oscillations 

kin 

in energy which can be explained due to breathing of 
the soliton after the collision. In the bottom panel E7- 
is given for two Vq. We find again that the velocity v t 
decreases for larger Vq. 
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FIG. 16: Density plot for V = 5.182 in a) and V = 5.183 in 
b) (hox = 40). In part a) the soliton decelerates and remains 
at the edges of the well before it continues to move to the 
right hand side. In b) the soliton slows down and remains at 
the edges of the well before it decides to move back to get 
trapped by the well. 
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Snapshots of the condensate density for different times 
are given in Figs [12] and [T3l The incoming soliton trans- 
mits almost completely through the well, only a small 
portion is reflected as radiation. Furthermore both fig- 
ures show once more that the radiation increases for 
deeper wells, i.e. the transmitted fraction is reduced. 

Next we look at the reflected part in Fig. [14J There the 
ratio between interaction and the density contribution to 
the kinetic energy is 



E L 



< i. 



(28) 



Comparing to Eq. (|2lj) this clearly indicates that the re- 
flected part in this regime is not soliton-like. Instead the 
almost vanishing absolute value for the interaction term 
shows that the main contribution, the kinetic energy, is 
being carried by radiation as proposed in Q. 

The findings in this section help to understand the fi- 
nite width of the transmission bands that we found in 
Sec. IIIII If Vi decreases the width of the transmission 
bands decrease as well, because less radiation is needed 
to trap the soliton. 

We illustrate the transmission and trapping behavior 
at the critical point for Vb in Figure \W[ 

In addition there is a second mechanism to trap a soli- 
ton that is similar to the temporary trapping of a linear 
wave packet, which occurs at the boundary between re- 
flection and transmission regions with radiation playing 
no role. For this situation we find that the soliton re- 
mains in the well for some time td until it reflects. We 
can measure this time delay td as the time the center of 
mass reaches the center of the well for the first time until 
it leaves the center again. Furthermore we see that by 
carefully adjusting the potential depth td can be large 
enough to observe a temporally trapped soliton in exper- 
iments. Fig. [17] shows the delay of the soliton during 
the transition through the well. It is remarkable that the 
whole soliton can be trapped with negligible losses due to 
radiation (> 99%). The losses are indeed much smaller 



FIG. 17: Time delay of the soliton due to the existence of the 
potential well (solid line). These are results for the simula- 
tions in Sec. IIIII Near the transition point from full reflection 
to full transmission (left vertical dashed line) the time delay 
increases very fast and becomes very large. This means the 
soliton remains in the vicinity of the trap on a very long time 
scale. Within the transmission region it decreases again until 
the trapping mechanism kicks in (right vertical dashed line). 
This picture for the nonlinear regime differs from the analyti- 
cally calculated time delay for the linear case (dashed-dotted 
line) not only in the position and value of its maximum but 
also for the nonexistent negative time delay (right of the dot- 
ted vertical line). 
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FIG. 18: Density plot for V = 4.842 (l box = 40). The soliton 
decelerates and more than 99% of the initial soliton is being 
trapped for a long time period (see also Fig. I17|l . This is 
connected to temporal trapping of a linear wave packet. 



than for the first trapping mechanism described before. 
The time evolution of the density given in Fig. [TH] shows 
an example of the temporal trapping of the soliton. This 
delay within the well is analogues to the interaction free 
case for a traveling wave packet with velocity v towards 
a well. There, an analytical expression for the time delay 
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is known [431 as 



d_ 

dE 



arctan 



1 



E 



2 y/E + Vo 



: tan 




E=v 2 /2 

(29) 



which is shown in Fig. [TT] as a dashed line. This delay 
can be explained as a temporary trapping of the linear 
wave packet during which it oscillates between both ends 
of the well before it escapes again. However, due to the 
nonlincarity the position and the value for the maximum 
time delay differs significantly from the linear case. Fur- 
thermore the time delay of Eq. ([29]) becomes negative, 
which happens if the quantum well is deep enough to 
turn a quasi-bound state into a bound state which is in 
contrast to the nonlinear case. One should however note 
that in the linear case about 42% of the wave packet re- 
flects at the point of maximal time delay. For smaller 
velocities this value seem to converge towards 50%. In 
addition the maximum time delay for the linear case lies 
well below the one for the nonlinear case, where it seems 
to diverge at the critical value for Vq. Therefore the con- 
nection to the nonlinear case is still unclear and needs 
further investigations that go beyond the scope of this 
work. 

We want to remark that although this would be an 
elegant way for lossless trapping of a soliton, Fig. [T7] also 
shows that the width of the this delay is very narrow and 
therefore harder to realize and to observe experimentally. 
In a BEC experiment with a small enough number of 
atoms it may be expected that superposition states will 
occur in this region [9j . Hence the other trapping method 
is favorable when it comes to experiments, even though 
one has to take into account minor losses. 



VI. PROBING ENERGY LEVELS 

Trapping of soliton amplitude is sensitive to bound 
states in the well. Data presented in the previous sec- 
tions has already suggested that trapping results from a 
resonant interaction of the soliton with a stationary de- 
fect mode, where the relevant energy scale is the soliton's 
chemical potential. By exploiting this resonant relation- 
ship, we are able to extract the bound-state energy by an- 
alyzing soliton scattering data. We proceed by compar- 
ing the scaled particle number of nonlinear bound state 
solutions with the trapped component after scattering a 
soliton with the same chemical potential. 

In this section we model the defect as an attractive 
delta potential V(x) = — Vq5(x), which has only one 
bound state at E b = -V r 2 /2. We solve Eq. © with 
g = — 1 with a soliton initial wave function ©, varying 
the amplitude A and thus the chemical potential fii = 
gA 2 /2. After the soliton has scattered, we integrate the 

scaled particle number Nl = f_ b \ip\ 2 dx of the trapped 
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FIG. 19: The scaled number of trapped atoms Nl vs. fii 
from time-dependent simulations of Eq. ([2]) is compared with 
Ns(fi) for stationary solutions from Eq. (130[) . Stationary so- 
lutions were also found with our numerical code for checking 
numerical accuracy. Results taken over a range of initial ve- 
locities Vi show a consistent picture. The most significant 
deviations occur at the onset of trapping around the location 
of the linear bound state at Eb = —0.5. This feature is most 
clearly distinguished for the smallest velocities. 



component (choosing b such as to capture at least 99% 
of the initial soliton's normalization). In Fig. [19] we com- 
pare this data with the normalization Ns = J \(j>\ 2 dx of a 
stationary localized solution ip(x,t) — 4>{x) exp(i/isi) of 
Eq. © with the same chemical potential fis = Mi- We 
analytically find 



N s = 2 V^2^ - V 



(30) 



which is shown as the full line in Fig. 1191 The energy 
E b of the linear bound state [of Eq. ^ with g — 0] is 
found at the intersection of the line with the fj, axis, i.e., 
E b = fi S at N s = 0. 

As expected, trapping is observed in the time- 
dependent simulation only for fa J$ E b (Figs. [TH] and 
|2T)]) with the scaled particle number increasing with de- 
creasing \ii , roughly following Eq. [30] As seen in Fig. [TO] 
where data with a variety of different initial velocities 
is compared, the trapped component is systematically 
about 20% larger than expected from the exact station- 
ary solution. We have verified that the final state of the 
trapped component in the time-dependent simulations 
corresponds to a stationary solution with further reduced 
chemical potential compared to the initial /i,. While at 
this time we are not able to explain why this happens, 
the important result is that recording the trapped com- 
ponent as a function of fa allows us to locate the bound- 
state energy E b . 

In Fig. [50] we have plotted the trapped component as 
a function of the soliton's initial chemical potential fa 
for different trapping potentials. Least square fits of the 
data (data points with Nl > 0.003 were included) to the 
functional form Nfn 



the bound state energy E b . 



fa provide estimates (3 for 
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FIG. 20: Scaled particle number of the trapped component 
Nl vs. after soliton-defect scattering as in Fig. [19] for dif- 
ferent values of the defect strength Vo . Short vertical lines in- 
dicate the energy of the bound state Eb for each of the values 
of Vb- Square-root fits to the data (as explained in the text) 
provide estimates (3 for Eb from the scattering data. Values 
for 13 found are -0.120, -0.505, -1.12, -1.97, -3.08, -4.43, which 
should be compared with the corresponding exact values of 
Eb given in the legend. 



We expect that bound state energy levels of narrow po- 
tential wells of more general shape than the one studied 
here could be probed experimentally by scattering bright 
solitons using this scheme. For defects with more than a 
single linear bound state, we expect that only the least 
strongly bound one can be detected in this manner. 



the defect size is small compared to or of the same order 
as the soliton length. We have found a rich transmission 
- reflection spectrum, which is strongly influenced by the 
level structure of the defect. In contrast to the scattering 
of linear waves, as in the scattering of single or indepen- 
dent atoms, part of the soliton can be trapped on the 
defect corresponding to the population of bound states 
or nonlinear localized modes. 

We find windows of transmission associated with 
above-well resonances in linear scattering. Nonlinear in- 
teractions modify the line shapes and lead to an abrupt 
onset of transmission. We have shown that particle loss 
of a few percent due to radiation leads to slowing down of 
the soliton due to energy conservation. A resonant cou- 
pling between the incoming soliton and bound states on 
the defect is identified as the mechanism leading to trap- 
ping and population of bound states. We have shown how 
this resonant coupling provides a way to experimentally 
probe bound states of an unknown localized potential 
well by scattering of nonlinear waves. 

In future work it will be interesting to study the sub- 
sequent scattering of two or more solitons or the scat- 
tering of solitons as we might expect additional effects 
of matter-wave enhancement or triggering the release of 
stored solitons. In addition, quantum many-particle ef- 
fects like condensate fragmentation or macroscopic quan- 
tum tunneling could be investigated in the context of 
resonant soliton-defect scattering. While this question is 
very interesting it leads beyond the scope of the present 
work and thus is deferred to further study. 



VII. CONCLUSIONS 

In this work we have investigated the scattering of a 
bright soliton on a linear defect in the context of matter- 
wave solitons. By numerical simulation and variational 
collective-coordinate studies, we have investigated the 
regime where the solitons are slow such that nonlinear 
energy scales dominate over kinetic energy and where 
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